Spin-current order in anisotropic triangular antiferromagnets 
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We analyze instabilities of the collinear up-up-down state of a two-dimensional quantum spin-5 spatially 
anisotropic triangular lattice antiferromagnet in a magnetic field. We find, within large-5 approximation, that 
near the end point of the plateau, the collinear state becomes unstable due to condensation of two-magnon 
bound pairs rather than single magnons. The two-magnon instability leads to a novel 2D vector chiral phase 
with alternating spin currents but no magnetic order in the direction transverse to the field. This phase breaks a 
discrete Z2 symmetry but preserves a continuous (7(1) one of rotations about the field axis. It possesses orbital 
antiferromagnetism and displays a magnetoelectric effect. 

PACS numbers: 



Introduction. The field of frustrated quantum magnetism 
has witnessed a remarkable revival of interest in recent years 
due to rapid progress in the fabrication and characterization 
of new materials and a multitude of theoretical ideas about 
competing orders and new quantum states of matter [ 1 1. Stud- 
ies of quantum 2D triangular antiferromagnets with spatially 
anisotropic exchange, such as CS2CUCI4 and Cs2CuBr4, are 
of particular interest because of their surprisingly rich phase 
diagrams in a magnetic field [|2][3]]. One of the best under- 
stood phases is a collinear up-up-down (UUD) phase with a 
fixed magnetization equal to 1/3 of the saturation value. In 
this phase, two spins in each triangle point up and one points 
down. The state breaks discrete Z 3 translational symmetry, 
but preserves continuous U(l) symmetry of rotations about 
the field direction. As a result, all spin excitations are gaped 

a. 

In a classical isotropic 2D Heisenberg systems with nearest 
exchange J, the UUD phase is the ground state for just one 
value of the external field h = 3 J ( 1/3 of the saturation field 
h sa ,t — 9 J). At all other fields spins order in a non-collinear 
fashion. In an anisotropic lattice with exchanges J and J' (see 
Fig. 1), a non-collinear order wins for all fields, so that clas- 
sically UUD phase is never a ground state. For quantum sys- 
tems, the situation is different as quantum fluctuations favor 
a collinear spin structure |5|. In the isotropic case, quantum 
fluctuations stabilize the UUD phase with gapped spin-wave 
excitations in a finite interval of h. In an anisotropic case, a 
large-S 1 analysis shows [6 1 that UUD phase persists up a finite 
anisotropy J' cr = J(l - ^3/105) (see Fig.[T}. The bound- 
aries of the UUD phase have been determined from the local 
stability analysis as the values of h at which spin-wave disper- 
sion softens. Of the two low-energy spin-wave branches, one 
softens at the lower boundary of the UUD phase and another 
at the upper boundary. Near the critical J'/ J, both spin-wave 
instabilities occur at finite momenta, and each leads to a chi- 
ral, non-coplanar state (distorted umbrella, in terminology of 
[6 1) in which (S r ) has finite components along both directions 
perpendicular to the field J6[ |7| . 

The analysis of the same model for S = 1/2, however, 
found very different states surrounding the UUD plateau near 




FIG. 1: (Color online) (a) Anisotropic triangular lattice with ex- 
changes J and J', (b) Distorted umbrella state, (c) Schematic phase 
diagram of the model in the vicinity of the UUD end-point at 6 = 4. 
Thin solid (red) lines mark single-particle instabilities of the UUD 
state at h c i, C 2(8). Thick solid (blue) line is the two-particle instabil- 
ity line towards a spin-current state, which emerges at 5 > <5 cr , and 
dotted (black) lines indicate phase transitions between the umbrella 
and the spin-current state. Dashed (red) line indicates a would-be 
single-particle instability, which is pre-empted by the two-particle 
instability. (Blue) arrows in the insert on the right show the arrange- 
ment of spin currents. 



its end point, which for 5 = 1/2 extends all way to J' = 
0. These states are collinear spin-density wave (SDW) 
states, with incommensurate spin modulations along the field 
direction but no long-range order in the transverse direction 
0. This discrepancy poses the question whether the phase 
diagram for S = 1 /2 is qualitatively different from the one at 
large S, or the ground states surrounding the UUD phase are 
different from the ones predicted by spin-wave theory even for 
large S. 

In this work we re-visit the large S analysis of the UUD 
state and show that the spin-wave phase diagram is incomplete 
for any S. We show that, prior to a single-magnon instability, 
the system undergoes a pairing instability, in which the two- 
particle collective mode, made of magnons from the two low- 
energy branches, softens at zero total momentum of the pair. 
As a result, the actual instability near the end point of UUD 
phase is towards the uni-axial state with no magnetic order 
in the transverse direction, similar to the situation for S — 
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1/2. We solve the "gap" equation for the two-magnon order 
parameter and show that it is purely imaginary. Such order 
parameter breaks a discrete Z 2 symmetry and gives rise to a 
bond-nematic state with non-zero vector and scalar chiralities 
within a single triangle of spins: (Sa ■ Sb x Sc) 7^ and 
(S A x S B ) = (S B x S c ) = (S c x Si) / (vector and 
scalar chiralities are proportional to each other since the total 
magnetization M — (S z ) is finite). Such a state supports 
circulating spin currents (Fig.[2]i and we label it a spin-current 
state (SC). We present the modified large-5 phase diagram of 
the model in Fig.[T] 

Experimental signatures of a SC state are rather peculiar. 
First, it exhibits a magneto-electric effect because both spin 
current and electric field are odd under spatial reflections and 
couple linearly iflOl . As a result, spin- wave excitations of the 
SC state depend linearly on E. Second, orbiting spin cur- 
rents generate charge currents, which in turn produce stag- 
gered magnetic moments, which can be measured by NMR 
and ^SRlfTTl. 

The model. We consider a system of localized spins 
on an anisotropic triangular lattice with Heisenberg nearest- 
neighbor interactions J and J', subject to an external field 

h = 2fi B H z : 



where at small k 



J'=l,2 



hSl 



(1) 



where ai, 2 = a(l/2, ±vo/2) connects spins on neighbor- 
ing chains, and a is the lattice constant. For convenience, we 
rescale h = hS and use h for the field. The saturation field, 
above which the magnetization M reaches maximum possi- 
ble value M sat = 5, is given by h s&t = (2 J + J') 2 / J. We 
are interested in the behavior of the system near /i sat /3, where 
quantum fluctuations win over classical fluctuations and stabi- 
lize UUD phase in a finite range of fields. In the isotropic case, 
,/' = J, the UUD phase exists in a field range between h c \ — 
(W3)(l-0.5/2S) andfr c2 = (/i sat /3)(l + 1.3/25). In the 
anisotropic case, J' < J, the competition between quantum 
and classical fluctuations is measured [6| by a dimensionless 
parameter 6 = (405/3)(l — J'/ J) 2 . As anisotropy increases 
and 5 grows, the width of the UUD state decreases and even- 
tually vanishes at S = 4, which defines J' — J' a . 

The excitation spectrum of the UUD phase at 5 < 4 can be 
straightforwardly obtained by using a three-sublattice repre- 
sentation for two spin-up and one spin-down sublattices and 
introducing ||6] |9) three sets of Holstein-Primakoff bosons, 
a, b, and c. One of the three spin-wave branches describes 
the precession of the total magnetization, has energy of the 
order h sat /3, and is irrelevant to our analysis. The other two 
branches describe low-energy excitations. Explicitly, 
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wi, 2 (k) = ±[h-h - —J- ~Jk 
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Z k = V 9 + 10S(6k 2 - 3(5fc2 + lOSk 4 ), 



Z k , (3) 



(4) 



and ho = J + 2 J'. The excitation di ^ softens at the lower 
boundary of the UUD phase, at h = h c i(S) — h en d — 
9 J/ (405)^(4- <5)/3, where h cnd = h Q (l + 17/(1205)). 
The softening happens at a finite momenta ±ki = (±/ci, 0), 
where hi « (3/(105)) 1 / 2 (l + x /(4- 5)/12). The exci- 
tation c?2,k softens at the upper boundary h — h C 2(S) = 
/i ond + 27 J/ (405)^/(4- S)/3, at momenta ±k 2 = (±fc 2 ,0), 
where fc 2 = (3/(105)) 1 / 2 (l - ^/(4 - 5)/ 12). The spin-wave 
softening at either h c i{8) or h C 2(S) signals condensation of 
one-magnon excitations. A Ginzburg-Landau-type analysis 
shows |6| that condensation spontaneously breaks Z 2 sym- 
metry between degenerate minima at ±ki and ±k 2 . As a 
result, one-magnon condensation gives rise to an incommen- 
surate spiral order with spontaneously broken 0(2) x Z 2 sym- 
metry and a finite non-coplanar long-range order (S^ ,y ) 7^ 0. 

At the end-point of the plateau 5 = 4, h cl = h c2 = h end , 
both spin-wave branches touch zero simultaneously at ±k = 
(±fc o ,0), where k = ^3/(105). The presence of four soft 
modes leads to a variety of possible non-coplanar chiral orders 
with non-zero (S^ v ). However, we show below that instead 
the system undergoes a pre-emptive pairing instability into a 
state with no transverse order, (S^ ,y ) = 0, but nonetheless 
with a finite chirality (z ■ S r x S r ') 7^ 0. 

Magnon pairing. To analyze a possibility of a bound state 
of two magnons, we need to include magnon-magnon interac- 
tion. The derivation of the interaction Hamiltonian is lengthy 
but straightforward: one has to express express two-magnon 
interaction Hamiltonian H„„ d , originally written in terms of 
Ok, & k and Ck bosons, in terms of the low-energy eigen-modes 
di } k and ci 2i k from Eq. Q. The full transformation is given 
in (9). Near momenta ±k , which are mostly relevant to the 
pairing problem, this transformation simplifies to 

a k =^(e is *d 1)k -e- iSk 4_ k ), 



/(k), 

c k = /(k)(d 2 ,k-e i2sk d t 1 .-k)- 



b k = -^( e -^d 1M + e^4_ k ), 



(5) 



(1 - <V4)fc 2 r 1/4 and 



where f(k) = y/ko~[{k x ± fc ) 2 
s k = 7r sign(fe a; )/4. 

Consider first 6 < 4, when only one boson becomes soft at 
either h c i or h c2 , while other remains massive and can be ne- 
glected. For concreteness, consider the vicinity of h c \, where 
d\ excitation softens. The magnon-magnon pairing interac- 
tion involving only d\ bosons is 



U 



(4) 
d\d\ 



8(J- 



(A-5) 



(6) 
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This interaction is obviously strongly repulsive and does not 
give rise to a bound state. The same holds for d,2 mode near 
h C 2- As a result, one-magnon condensations at h c i and h C 2 are 
the true instabilities, and the system develops a non-coplanar 
spiral order at h > h C 2 and h < h c \. 

For S « 4, the situation is different. Magnon-magnon inter- 
actions within d\ or d,2 sectors are still repulsive, but now we 
also have interaction between d\ and e?2 bosons, both of which 
are gapless at ±k . The di — d 2 interaction with zero total 
momentum has two relevant terms: one describes "normal" 
2—^2 process with simultaneous creation and annihilation of 
d\ and bosons, the other describes "anomalous" 4 — > and 
— > 4 processes with simultaneous creation or annihilation of 
two di and two di bosons. We find that the strongest pairing 
interaction involves momentum transfer ±2fc for each of the 
bosons involved. The corresponding interaction reads 

n dld 2 = ^H^^'^^Uo+P^l-ko-p^l.-ko+q^.ko-q 

-4 v 4-„dl w „4 w , „4 w + h - c - ( 7 ) 

l,k +p 2, — k — p 1, — k +q 2,k. — qJ v ' 

where p and g are much smaller than k , and the vertex 

*(p,«) = -(J + 2J')/ 2 (p)/ 2 (9) -(J + 2j')|i (8) 

where /(p) was introduced after Eq. ((5), and the limit 
stands for 5—^4. The pairing interaction with small mo- 
mentum transfer, g)4,k + P 4 -k Q - P d i,k +q^2 ,-k -q, 
has a much smaller $(p, g) which remains finite in the limit 
p, q — >• 0. Such interaction is then irrelevant for our analysis. 

Now observe that the sign of 2 — s- 2 term is negative, while 
the one of 4 -> term is positive. The negative sign of the 
2^2 term implies that the "normal" interaction between d\ 
and di bosons is attractive and favors a pairing with 

^k (p) =< di, ko+ prf 2 ,-k -p >= T/|p| = F—k a ( P ). (9) 

The negative sign of the 4 — > term does not allow the so- 
lution with real T (the corresponding coupling constant van- 
ishes), but instead favors a solution with imaginary T = iT. 
For such T, the pairing vertex which couples to 4 — » term 
has opposite sign compared to the vertex which couples to 
2—^2 term, and this extra sign change compensates the sign 
difference between 2 4 2 and 4 — !• interactions. Note that 
since the Hamiltonian |7]) does not conserve the number of 
bosons, the order parameter does not possess a £7(1) phase 
symmetry. In practice, this implies that the gap equations for 
real and imaginary T's are different. And, in fact, the symme- 
try that is spontaneously broken at the transition is Z2, corre- 
sponding to the sign of T. 

For T = iT, the linearized "gap" equation reads at 6 = 4, 

NS^ p 2 w 1 (k + p) + w 2 (k + p)' 

(10) 
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FIG. 2: (Color online) The structure of spin currents in the SC state. 
The domain wall, denoted by vertical (red) dotted line, separates do- 
mains with opposite chirality. 

Substituting the dispersions, we find 

i = IAyA (id 

p 1 1 

It is important that the integrand scales as l/|p| 3 , so that the 
2D integral over p diverges and overcomes the smallness of 
1/S in the pre-factor. In l/|p| 3 , one power of l/|p| comes 
from the dispersion and the other two powers are due to the 
divergence of the coherence factor f(p) at p — > 0. Away from 
(5 = 4, |p| is replaced by (|p| 2 + (1 - <S/4)/s 2 ) 1/2 , and the 
integral in the r.h.s of dllfc behaves as l/\/4 — S. Collecting 
powers of 1/S, we find that a nonzero T emerges at i5 CI = 
4-0 (1/S 2 ). 

For completeness, we also analyzed possible pairing with 
the total momentum ±2k , but found that there is no enhance- 
ment of the kernel of the gap equation by coherence factors 
and, hence, no instability at large S. 

Spin-current order. The two-magnon instability does not 
lead to a conventional spin order in the direction perpendic- 
ular to the field because (di t k) = {d2,k) = 0. Fk (p) ~ T 
does not give lead to modulations of or the bond order be- 
cause the condensate does not contribute to magnon density 
or to (Sa ■ Sb) |9|. However, one can easily verify that for 
each triangle we now have (z ■ Sa X Sc) = (z ■ Sc x Sb) = 
{z ■ Sb x Sa) oc T, which implies a finite vector chirality 
and orbital spin currents which run in opposite directions in 
neighboring triangles, Figure [2] Note that the sign of Ising 
order parameter T determines the sense of spin current circu- 
lation. In our case vector chirality generates a non-zero scalar 
chirality (S^ - Sg x Sc) ~ T as well, because of the finite 
magnetization M along the z (magnetic field) axis. For tri- 
angles separated by distance r, z ■ (S(0) x S(r)) scales as 
Tcos k r. 

A spin-current (SC) order in dimensions D > 1 is nor- 
mally associated with non-coplanar spin ordering when the 
spins spontaneously select the direction of rotation in the XY 
plane. Remarkably, in our case the SC order appears in the 
absence of the standard spin order in the XY plane. 

The emergence of the SC order can be thought of as spon- 
taneous generation of Dzyaloshinskii-Moria (DM) interac- 
tion. Indeed, the interaction Hamiltonian (|7]i can be written 
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as HZ 2 = -(VJ/N)H™ H™, where 

r 

= i E fk (di,kda,-k - <k4,-k) ■ < 12 > 
ke±k 

As a result, the development of a non-zero T can be 
viewed as the appearance of Dzyaloshinskii-Moria interaction 
D(H° M +H D f o ), with D ~ T. This observation helps to un- 
derstand magneto-electric effect in the SC state: because D is 
a pseudoscalar, it couples linearly to an electric field E, i.e., 
D = Dq + D\E + .... As a result, spin-wave excitations of 
the SC phase depend linearly on E. 

SC order has been previously explored in ID spin ladders 
lfT2UT4 l and was suggested for a frustrated Heisenberg model 
in 2D [15]. There, however, a SC state is a spiral state, 
in which a continuous £7(1) symmetry is restored by strong 
quantum fluctuations. In our case spiral states are present in 
the phase diagram away from the end-point of the UUD phase, 
while the SC state emerges as a result of a pre-emptive two- 
magnon instability rather than due to divergent one-magnon 
fluctuations. Our two-magnon instability (which necessary 
leads to an imaginary order parameter) is also fundamentally 
different from two-magnon instabilities with real order param- 
eter which lead to a spin-nematic order, either on a site or on 
a bond |[T6lfT9l . Such order generally occurs in systems with 
ferromagnetic exchanges at least on some of the bonds, when 
there is an attractive interaction between magnons. Here, 
all exchange couplings are antiferromagnetic, and magnon- 
magnon interaction is repulsive. Our pairing of magnons from 
different branches is conceptually similar to the inter-pocket 
pairing in multi-band fermionic systems, such as Fe-based su- 
perconductors with only electron pockets [20]. 

The phase diagram near the end point of UUD state has 
been recently analyzed in |7| in a self-consistent semiclassi- 
cal formalism. This method, however, does not allow for the 
analysis of two-particle instabilities. 

Comparison with SDW state. Although our analysis 
uses 1/5 expansion, it is nevertheless instructive to compare 
symmetry properties of our spin-current state with that of a 
collinear SDW state observed for 5 = 1/2 near the end point 
of the UUD phase. Like we said, spin-current state is much 
closer to SDW state than a spiral state (the result of one- 
magnon condensation) because both spin-current and SDW 
states preserve U(l) symmetry of rotations about the field di- 
rection. But the two states do differ as SDW state has no 
chiral order [8|. It may be that 5 = 1/2 is simply spe- 
cial and non-chiral SDW state is only present at 5 = 1/2. 
But it also may be that the two-magnon instability, which we 
found, is only a 'tip of the iceberg', and the two-magnon con- 
densation triggers the development of multi-magnon conden- 
sates at some 8 > <5 cr , which in turn changes the properties 
of the spin-current state. This last possibility is inspired by 
the observation that SDW state is incommensurate and that 
the UUD-SDW transition for 5 = 1/2 is a commensurate- 



incommensurate transition 0. Such transition occurs via a 
proliferation of solitons - strings of displaced spins which are 
shifted from their equilibrium UUD pattern. Since chang- 
ing the direction of a single spin 5 requires 25 magnons, a 
proliferation of solitons implies condensation of 25 magnons 
per every displaced spin. Then, in magnon description, a 
commensurate-incommensurate transition involves a conden- 
sation of an infinite number of magnons. One can imagine, by 
analogy with coupled superconducting and spin density orders 
ll22l . that proliferation of SC domain walls, depicted in Fig. [2] 
may cause the appearance of an incommensurate modulation 
of (5 Z ) due to "density-density" type coupling between the 
magnon density and the density of domain walls. Whether or 
not this is the case requires going beyond the instability con- 
dition ( fTTj ) and analyzing excitation spectrum and inter-pair 
interactions within the spin-current phase 1 23 1 . 

Conclusions. We have described a novel two-magnon 
pairing instability of the up-up-down phase of the spatially 
anisotropic triangular lattice antiferromagnet in a magnetic 
field. The magnon pairing is of "inter-band" type in that the 
condensate is made out of bosons from the two different spin- 
wave branches. This instability pre-empts a single-magnon 
condensation for arbitrary spin 5 and gives rise to a highly 
unconventional 2D order in which transverse spin components 
are disordered, yet the ground state has a non-zero vector chi- 
rality on every lattice bond and circulating spin currents in ev- 
ery elementary triangle. This state breaks Z2 chiral symmetry 
but preserves U(l) symmetry of rotations about the field di- 
rection. The development of such a phase can be thought of 
as a spontaneous generation of the Dzyaloshinskii-Moriya in- 
teraction. 
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SUPPLEMENTARY MATERIAL FOR "SPIN-CURRENT 
ORDER IN ANISOTROPIC TRIANGULAR 
ANTIFERROMAGNETS" BY A. V. CHUBUKOV AND O. A. 
STARYKH 

One-magnon excitations in the UUD phase 

One-magnon excitations in the UUD phase in the 
anisotropic case (J' < J) have been analyzed in Ref. Q] For 
completeness, we present here the details of the derivation. 
We will use some of intermediate formulas in the next section, 
when we derive the pairing interaction between magnons. 

Spin-wave description of the UUD state proceeds as fol- 
lows. First, we use a three-sublattice representation where 
spins point up on sublattices A and B and down on sublat- 
tice C, and introduce Holstein-Primakoff bosons a, 6, and c 
respectively. Spins on the A sublattice are described by 



S Z A (V) = 5 - dOr, 

S+(r) = V2S\l- 



ala T 



^(l 



a].a r 

~4S 



K,(l) 



and spins on the B sublattice are represented by the same ex- 
pressions with a r replaced by 6 r . The expansion of the square- 
root is valid for 5^1, and below we assume that 5 is indeed 
large. Spins on C sublattice points opposite to those on A and 
B sublattices, and we have 



S* c (r) = S + 4c r , 5 c (r) = V2S(1 - ^)c r . (2) 

Plugging this in ([TJ, we obtain spin-wave Hamiltonian as the 
sum of the linear (harmonic) term and the interaction 



terms 

h{2) = S Y. [h^k + 7k(6 k 4 k + ct a t_ k ) + h .c] 

k 

+/i(a k a k + 4&k) + (2h - h)c{c k . (3) 

Here the sum extends over the magnetic Brillouin zone, whose 
area is 1/3 of the total area of the Brillouin zone, 7k = 

Je lk * + 23' cos(V3k y /2)e- lk */ 2 and h = J + 2J'. The 



interaction term ffW = JjW + H ™ is the sum of the trans 
verse (_L) and longitudinal (z) contributions 



H 



(4) 



("j) 77 51 {7i(444 a i+2+3 +&I44 c i+2+3 

ki-k 3 



+446361+2-3) + l-i(b[ ala 3 a 1+2 -3 + <444 c i+2+3 



!,t„t 



44&3&1+2+3) +h.c.}, 



(4) 



and 



77 {7l-24 a 26^6i_ 2 +3 



-7l_ 2 6[6 2 4 C l-2+3 - 72-l4 fl 24 C l-2+3}- (5) 

Here we denote for brevity 71 = 7 kl , a\ = and so forth. 

If we take only the quadratic part ^ and diagonalize it, 
we find that the UUD phase is stable for just one field value 
h = 3 J for the isotropic case (J = J') and is unstable for all 
fields when 5 7^ 0. (In the notations which we used in the text, 
the degree of anisotropy is measure in terms of dimension- 
less parameter 6 = (405/3) (1 - J'/J) 2 , Ref 1). However, 
interactions between magnons stabilize the UUD phase over 
the finite S range. To see this, we first modify the quadratic 
form by adding the leading 1/5 Hartree-type self-energy cor- 
rections from 4-boson interaction terms Q and (BJ, and then 
diagonalize the effective quadratic Hamiltonian. Since we are 
interested in large 5 and small anisotropics (when S = 0(1)) 
and in field range near h = 3 J, Hartree corrections can be 
computed at S = and h = 3 J, when classical UUD state is 
critical and its spin-wave excitation spectrum does not contain 
complex modes. In 2D Hartree corrections are all finite, and 
adding them to <|3j we obtain the harmonic Hamiltonian of the 
UUD state in the form 

-ffuud = [[7i,k4 & k+72,k(44 k +44 k )+ h -c-] 



-(h + £i)(<a k + 6 k 6 k ) + (2h + S 2 - /i)4c k .(6) 



where 



7j,k = 7k + 



(7) 



and the self-energy components are Si = 0. 14 J/5, S 2 = 
0.67 J/5, Si j0 = -0.11 J/5 and Yl 2 = 0.18J/5. Observe 
that |6]) reduces to Q in the 5 —> 00 limit. 

Alternatively, one could first diagonalize the linear spin- 
wave part Q, express the 4-boson interaction part in terms 



2 



of new operators of system's eigen-modes and then correct 
magnon dispersion by adding to it quadratic terms in new op- 
erators, which appear as a result of normal-ordering of Eqs. 
Q and Q in the new basis. Such a procedure was first ap- 
plied by Oguchi and is known as Oguchi's corrections |2] 

Low-energy excitations near k = encode the important 
physics, and in this region analysis of Eq. |6]) simplifies con- 
siderably. Here we have 7^ ps 7 j ^ + ir k , with 



7j,k = h + o - -.Jk 2 , T k = ( J - J')k x . 



(8) 



Diagonalization of Eq. (|6]l proceeds in two steps. We first 
diagonalize H' mA , which is obtained from |6} by setting I\ = 
0. This is done by introducing new operators 



Pk = (flfe - h)/V2, d k = (0* + b k )/ \f2 



(9) 



which decouple in H^ ud . However, d-mode couples to c- 
boson via % k(dtc_t + h.c.) term, and to diagonalize the full 
quadratic Hamiltonian one needs to apply rotation 



dk = cosh 9 k u k + s'mh 8 k v^_ k 



where 

tanh(26» k ) = 



sinh OkUk + cosh O^v 



-fcj 



-2V27 



2.k 



2y/2 



(10) 



(11) 



2h + Si + £ 2 + 7i,k 3 

The quadratic Hamiltonian in terms of Uk, Vk andpk operators 
is 

H Ld = s y^[^ P P k Pk + ^i,w k «k + w„tt k u k ]. (12) 

k 

The u boson describes the precession of the total magnetiza- 
tion. The corresponding frequency ui u (0) ~ h is large, which 
implies that this mode is irrelevant for low-energy physics. 
The two remaining bosons, p and v, are the low-energy modes 
of interest. For small k 

w„(k) = [h-h° cl +2(J- J')] + ^Jk 2 , 

w„(k) « [/i° 2 -2(J- j')-/i] + ^Jfc 2 . (13) 

where = 3J - 0.5 J/ (25) and h° c2 = 3 J + 1.3 J/(25) are 
the boundaries of the UUD phase in the isotropic case. 

The next step is to account for the remaining part of i? uu( j 
in ([6j, which is proportional to T k = ( J — J')k x . The corre- 
sponding term, which we denote by H" ud , has the form 



^u'ud = 3^^r k ( Pk «_ k -h.c), 



(14) 



The diagonalization of H' uud + H" d proceeds in the same 
way as before: we introduce new operators e?i. k and d2,k as 



Pk = cosh 4>kd\.k + i sinh 4>kd 2 _ k , 
v_ k = cosh k c?2._ k + i sinh <f>kd\ k , 



(15) 



and choose (pk to eliminate non-diagonal G?i. k d2.-k terms. 
This last requirements leads to 



tanh(20 k ) 



6{J-J')k x 



6(J - J')k x 



u p {\s)+u v Qs) Ah + 3 J{k% + h*) 



(16) 

Here Ah = h° c2 - h° cl = 1.87/(25) is the width of the UUD 
phase at J' = J. The diagonalized quadratic Hamiltonian is, 
up to a constant, 

-Huud = S y^i^i.k^-.k + U2dl M d 2 ,k], (17) 

k 

where at small k 

wi, 2 (k) - ±±(lj p - u v ) + \^{u p + u v ) 2 - 36(J - J') 2 k x 



= ±[h-h Q -^J-^jtf 



205 



(18) 



with Z k = ^9 + 105(6k 2 - 35fc 2 + 105k 4 ). 

The UUD phase is stable with respect to small perturbations 
when both modes are positive. The full analysis has been done 
in Ref. 1 where it was shown that UUD phase survives up to 
<5 = 4. For our purposes, we focus on the region near the end 
point. The UUD phase is stable at h C 2 > h > h c \, where 



h c x = h 



end 



9J 
405 



-,h c2 = h 



end 



27J 
40S 



(19) 

and /i cn d = ho(l + 17/1205). Near the lower critical 
field h c i, the mode u)i(k) softens at ±ki = (fci,0), where 
ki « (3/(10S')) 1 / 2 (l + ^(A- S)/12). Near the upper criti- 
cal field h C 2, the mode ui2(k) softens at ±k2 = (k 2 , 0), where 
fc 2 « (3/(105)) 1 / 2 (l- v /(4 - 6)/ 12). At 5 = 4, the two crit- 
ical fields become equal h c i = h C 2 = /i en d> and both modes 
soften at the same k\ = k 2 = fc 2 = 3/(105). At this 5, the 
excitation spectra are 



Wl(fc) 

w 2 (fc) 



3J 

~2 
3J 

T 



(A'.? 



fc 2 ) 2 



4fc 2 fc 2 



(k 2 . 



3J 

T 
3 J 



{kl 



fc 2 ) 2 +4fc 2 fc 2 + T (fc2_fc 2 ^20) 



Observe that at this point (J - J')k Q = 187/(105") = 2 Ah 
and 3Jfcp = Ah, hence tanh(20/ £Oi o) = 1 in (16 1, i.e., k 
diverges at ±ko = (ko,0). The divergence of <fi implies that 
the coherence factors cosh^fc and sinh</>fc strongly diverge 
too. At small deviations from ±kp and 5 = 4, 



cosh 20 k ss sinh 20 k = 
2fcp 



1 



= / 2 W (21) 



(fc 2 -fc 2 ) 2 +4A : 2 fc 2 + (4-5)fc 4 



Below we will need to express bosons a, b, c via the low- 
energy eigen-modes d\ and d 2 . Working backward through 



3 



transformations Q, (jTUJ, and (15) , we obtain 
1 



6k = 



72 



{(cosh 0k + isinh 0k)di,k 



— (cosh 0k — isinh^k)^ -kl> 

-1 



\/2 



{(cosh 0k — « sinh </>k)rfi,k 



+ (cosh 0k + « sinh 0k)<4,-ki' 
Ck = -\/2{cosh0kd2,k — isinh^k^i _ki- (22) 



Near k ± ko and i5 = 4, 0k is large, and using (21 1 one can 
simplify the transformation to 



V2 



(e ZSk d 



l.k 



4,-k), 



/(k) 



(e" lSk di,k + e 4Sk 4_ k ), 



c k = /(k)(rf 2 ,k-e i2s M t 



1,-k 



)• 



(23) 



Here Sk = 7rsign(fc x )/4. 



Derivation of the pairing interaction between di and di 
magnons 



can be expressed as 

3i7^ 3 / 

U ( Z 2 = £ ((! + cosh2 ^ cosh2 ^ + 

+ sinh 20 9 sinh 20 p ) x (^ R (q)V L (p) + V{(q)* R (p)) 
+(1 + cosh 20 g cosh 20 p — sinh 20 9 sinh 20 p ) 
x(* R (q)* R (p) + *\ j (q)* L (p)) + 
+(—1 + cosh 20 g cosh 20 p + sinh20 ? sinh 20 p ) x 
x(* R (q)*{(p) +h.c.) 

+(—1 + cosh20 g cosh20 p — sinh20 9 sinh20 p ) x 
x(y R (q)y R (p) + ^ L (q)^ L (p)+h.c.)). (24) 



We see that there are two types of pairing vertices: the ones 
with transferred momentum (for a given boson kind) of the 
order 2fco (these are ^l'&r. terms), and the ones with trans- 
ferred momentum near zero (^l^l and ^r^b terms). For 
the first set of terms, the vertex contains cosh 20 g cosh 20 p + 
sinh 20 g sinh 20 p and diverges at 5 = 4 in the limit p, q > 0. 
For the second set, the vertex contains cosh 20^ cosh 20 p — 
sinh 20 g sinh 20 p , and the leading divergent terms cancel out. 
As a result, the vertex with momentum transfer near 2ko is 
much stronger. Keeping only this vertex and using the asymp- 
totic forms of cosh20 g / p and sinh20 g / p from Eq. (21 1 we 
obtain 

3J 3 



H 



(4) 
d 1 d 2 



2 

q :P 



r(?)(**(q)-*k(q)) 



x/ 2 (p)0Mp)-* f L (p)). 



(25) 



To obtain the interaction between low-energy magnons, one 



has to express the interaction Hamiltonian = H, 



(4) 



written in terms of a, 6, and c bosons, Eqs. |4} and ((5), 



via di and c?2 operators with the help of Eq. ( |22| i, and find 
which of the generated interaction terms are the strongest. 
This procedure is straightforward but time-consuming. We 
analyzed pairing interaction with zero total momentum of the 
pair and with total momentum 2fco- We found that the inter- 
action matrix elements are much stronger for the former case 
(zero total momentum pairs). The computational procedure 
is similar in both cases and we present only the details of the 
derivation of the strongest interaction. 

Because fco = (3/(10S')) 1 / 2 is small, we approximate the 
factors 7fe in Eqs. Q and <|3j by their values at k = 0, i.e., 
approximate 7fc by 70 = J + 2 J' ps 3 J. We verified that that 
keeping the momentum dependence of 7k's only gives rise to 
irrelevant small corrections. 

We assume and then verify that the dominant contribu- 
tion to magnon pairing comes from momenta near ±ko. 
To obtain the pairing vertices with zero total momentum, 
it is then convenient to introduce pair operators ^/{(q) = 



hq^2,- 



and * L (q) 



*1 -ko 



hq d 2,k -«< where 



|q| << fco. Expressing iJ^ 4 ' in terms of d\ and e?2 we find af- 
ter long but straightforward calculation that the pairing vertex 



Solution of the gap equation 

As is customary in superconductivity studies, we add to 
the Hamiltonian infinitesimally small pairing terms Ql = 
^o,L(q)^L{q),QR = $o : i?.(<?)* h{q) with generally complex 
$o,l an d &o.R an d obtain the renormalized $l(<7) and $#(9) 
by summing up ladder series of vertex corrections. At the 
pairing instability, the pairing susceptibility diverges, and the 
equations for $l( < ?) an d ^ii(<z) have solutions even when we 
set bare $o.l(<7) an d $o.r(q) to zero. 

The diagrams for the fully renormalized <f>L(<z) and &R(q) 
at the instability are shown in Fig.[T] One can easily make sure 
that the full set of coupled equations for different $ and <£>* 
separates into two independent sets for $x and $^ and for $^ 
and &* L . Because the momentum dependence of the pairing 



vertex in Eq. (25 1 is factorized into f 2 {p)f 2 (q), we search 
for the solution in the form $l(<7) = T/ 2 (q). Substituting 
this form into the diagrams and using one-magnon dispersions 
from ( |20] > we obtain after a simple algebra 

3J/ 4 (P) 



T 



(T - T*) 



NS ^ wi(fc 



-p) + cj 2 (fco +p) 

k 

{p 2 + (1 - 8/4)k$) 3 / 2 



(26) 



4 



/ = 



A 




FIG. 1: (Color online) Coupled set of diagrams for anomalous ver- 
tices $ l (upper line) and <£>|{ (lower line). 



We used the fact that (fcn + p) +uj2(k n +p ) = 6Jk 2 /f 2 (p) 
(see Eq. (20 1). It is obvious from Eq. ( 26 1 that T should be 
purely imaginary, T = iT. Substituting this into (26 1, we find 



that the equation for T has a non-trivial solution when 



" ~ "2^ |p|3 

p 1 1 



SN 



Generalizing this to 5 < 4 case, and replacing the sum over p 
by integral we find that the condition on the pairing instability 
reduces to 



(fp _ _S_ 

(p 2 + (1 - <5/4)fc2)3/2 



(28) 



where a = 0(1). Evaluating the integral we find that the two- 
magnon instability occurs at 4 — S = 0(1/ S 2 ). 



tions of motion for pair operators. We obtain 



id t v R (k) = n k * fl (k)--53*(p,k)(*t(p)-* i ( P )) 

p 

id t * L (k) = n u ^(k)-^^$(p,k)(^(p)-* R (p)) 

p 



n^ R (k) *(p- k)(*i( P ) - * i(p )), 

p 

n k *t.(k) _ A £$(p,k)(*t fl ( p ) - * fl ( P )), 
p 

(31) 



where f2 k = 5(wi(k + q) + W2(-ko - q)) and $(q,p) ~ 
—3Jf 2 (p)f 2 (q). We now Fourier transform i-dependence ( 
id t —> lu) and set to — because we are seeking the condition 
for the pair condensation. We then take expectation values of 
both sides of equations and form appropriate linear combina- 
tions to obtain 



(27) (v^(k) - * H (k)) 



(*t (k) - * L (k)) 



6J/ 2 (k) 3 
ft k AT 

6J/ 2 (k) 3 



p 

E/ 2 (p)(*Up)-^(p))- 

p 

(32) 



where we used 7^ w 3 J. We immediately see that and 
$l must be purely imaginary, and $>L,n(k) oc f 2 (k). The 
self-consistency condition then gives 



3 v 6J/ 4 (p) 1 3_ v _fco_ 



N 



(33) 



which is the same condition as Eq. (|27]>. 



Alternative derivation ofEq. (|27[). 



Pair condensation as a result of spontaneous generation of 
Dzyaloshinskii-Moria interaction 



We start from simple observation that 

[ttffi,*K(q)] = -S(wi(ko +q) +^ 2 (-k - q))*K(q), 

(29) 

and similarly for ^(q). We also see that 

[*fl(q),*fl(p)] = ^, J >(l+4,k + q ^l,ko+q+4,k - q ^,ko-q) 

(30) 

and [^^(q), ^(p)] = because for q,p <C ko, the two pairs 
do not overlap in momentum space. Inside the UUD plateau 
(<iL 2 k <ii/2,k) = 0, where the average is over the ground 



state. Hence the right-hand side of ( |30| l can be replaced by 
5 q .p, implying canonical bosonic commutation relations for 
pairs ^r/l- This allows for an easy derivation of the equa- 



We first observe that Eq. (|2"5j) can be re-written as 



(4) 



where 



< M = *E/ 2 (q) (**(q)-*k(q) 

q 

^™=*E/ 2 (q) (^(q)-^(q) 



(34) 



(35) 



Note that the integrands in ( 35 1 is the same as the two-magnon 
order parameter 



* J2 / 2 (q) ( rf i,q d 2,-q-dL q 4 



2,-q 



(36) 



qS±ko 



5 



Hence, once two-magnon condensation occurs, (H m 



±fco< 



ac- The condensate emerges at q = — k, and we obtain 



quires a non-zero expectation value, proportional to T, i.e., 
the Hamiltonian acquires an extra term 



n 



(4) 
d 1 d 2 



D [H 



DM 

+k 



u 



DM" 

— No- 



where 
D 



(37) 



(38) 



(See next section for a very similar calculation.) We now com- 
pare the Hamiltonian in Eq. ( |37| ) with the one which describes 
Dzyaloshinskii-Moria (DM) interaction in a triangular magnet 
and show that they are identical. The DM interaction on a tri- 
angular lattice reads [ 3 ] 



(Crflr) = -^^/ 2 (k) e - i ^<d 1 ,_ k d 2 , U -4 i _ k 4 !k ) 



^/ 2 (k)e— (zT-(*T)*)/ 2 (k) 



V2N 



2\/2zTcos^l4XN 



// 2 



4 J iV ^ k 2 x + k 2 + (\- J/4)fcg 



(44) 



Ht>m— 



DM-latt 



z ■ (ScW x S B (r) + S A (r) x S c (r) + 

r 

+S B (r) x S c (r)Y (39) 



In our case z coincides with the direction of external field. 
Expressing the spins in terms of Holstein-Primakoff bosons 
a, b, and c and transforming to d\ and d 2 operators using ( |23] > 
we obtain that after some algebra d22} gives 



Herefc is the deviation from ±ko, near which the phase factor 
takes values Sk = ±7r/4 correspondingly. Because (c r a r ) is 
imaginary, it does not contribute to S^-S^ cx Re (S 1 ^ 5^), but 
the spin current Jq A becomes non-zero. Similar calculation 
for other bonds shows that (c T a T ) = — (c r b r ) = (aJ6 r ). These 
relations fix the relative signs of spin currents and lead to two 
current patterns shown in Figure |2]of the main text. 



DM-latt 



6S(H 



DM 
+feo 



n 



DM • 

— fcn , 



(40) 



Comparing with ( 37 1 we immediately see that the appearance 
of a two-magnon condensate with an imaginary amplitude T 
can be viewed as a spontaneous generation of DM interaction 
with the coupling D cx T. 



Structure of spin currents 

The z-component of the spin current on the bond (n,m), 
connecting sites n and m, is defined as 



72 _ _ / c- C+ _ C+ C - 



(41) 



Re-expressing the r.h.s. of this expression in terms of a, b, 
and c bosons, we find that spin currents along the bonds be- 
tween spins from A, B, and C sublattices belonging to the 
same elementary triangle at a coordinate r are determined by 
the following combinations 



JcA ~ iS(c^a£ — c T a r ), Jqb 
J 7 Xb = iS{b+a r - a+6 r ). 



iS(c^b^ 



C r b r ), 



(42) 



It is easy to generalize this calculation for the spins located 
at distance R apart from each other (we assume that R 3> 



(z-Sc(R) x S A (0)) cxIm(cRao) 

(•OO /»27T 

dk 



f 

Tkl cos[k • R] / 
Jo 

Tfcg cos[ko • R] / 
Jo 



Tkl cos[k • R]X (i?/0, 



ke 



ikR gos[' 







(l-5/A)k 2 



dk 



kJ {kR) 



k 2 + (1 - 5 /A) k 2 



(45) 



where Jo and Kq are Bessel functions and £ _1 = 
ko \J\ — 5 /A. The correlation decays exponentially for R ^> 
£. This implies that £ has the meaning of the radius of a 
bound state. Using the relation i5 cr = 4 — 0(1/ S 2 ), we obtain 
^S/k . 



Using (|23]l, we obtain 



(c r a r = 



x(e^di, q 



•(k+q) 



/(k)/(q)((d 2 ,„ 



„i2s k jt 



"1,-kJ 



(43) 
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